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SUMMARY 

>  Experimental  performance  and  blade  stresses  measured  on 
a  three-bladed  VTOL-type  propeller  tested  in  free  air  are  presented. 
The  isolated  propeller  was  tested  over  ranges  of  prop  speed,  forward 
velocity,  blade  angle  setting,  and  thrust  axis  to  free-stream  angle. 

Correlation  with  a  theoretical  method  of  predicting  propeller 
performance  and  blade  stresses  was  made.  It  was  found  that  when  the 
propeller  was  operating  in  a  flight  condition  for  which  the  theory  was 
developed,  correlation  between  theory  and  experiment  was  good.  When 
the  propeller  was  operating  in  a  flight  condition  where  very  small  posi¬ 
tive  or  negative  effective  angles  of  attack  were  encountered,  correlation 
between  theory  and  experiment  was  poor. 
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FOREWORD 

The  technical  effort  reported  herein  was  conducted  at  Cornell 
Aeronautical  Laboratory,  Inc.,  during  the  period  from  June  1965  to 
January  1967.  Mr.  Andrew  R.  Trenka  was  the  project  engineer. 

This  program  was  sponsored  by  the  U.  S.  Army  Aviation  Materiel 
Laboratories,  Fort  Eustis,  Virginia,  and  was  administered  by 
Mr.  John  E.  Yeates. 

Recognition  is  given  to  Dr.  I.  C.  Statler  and  Mr.  R.  P.  White,  Jr., 
for  their  assistance. 
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INTRODUCTION 


Recent  efforts  to  develop  V/STOL  aircraft  prototypes  have 
demonstrated  that  existing  means  for  predicting  propeller  performance 
in  the  static  and  transitional  flight  modes  are  inadequate.  Attempts 
to  alleviate  this  situation  are  being  made  under  U.  S.  Army  spon¬ 
sorship  at  Cornell  Aeronautical  Laboratory,  Inc.,  (CAL)  (Refer¬ 
ence  1)  and  at  other  research  establishments  (Reference  2),  These 
recent  advances  characteristically  involve  the  use  of  high-speed 
digital  computing  machines  to  carry  out  the  necessary  calculations. 

The  method  for  predicting  propeller  performances  and 
blade  stresses  developed  at  CAL  (Reference  1)  can  treat  the  static, 
transitional,  and  axial  flight  cases.  Included  in  this  mathematical 
model  are  the  effects  of  wake  contraction  (for  the  hovering  case), 
consideration  of  finite  number  of  blades,  provisions  for  incorporation 
of  experimental  sectional  airfoil  characteristics,  effects  of  structural 
deformations,  and  first-order  effects  of  wing  and  nacelle  interfer¬ 
ences  on  the  blade  aerodynamic  loading. 

During  a  previous  program,  an  attempt  was  made  to  correlate 
the  CAL  theory  with  data  collected  at  NASA  Ames  in  the  40 -foot  by 
80-foot  low-speed  wind  tunnel.  Data  which  were  obtained  on  a  0.  6 
scale  model  of  the  XC-142  are  presented  in  Figure  18  of  Reference  1. 
The  experimental  error  associated  with  the  measurement  of  the  thrust 
and  torque  during  those  tests  was  estimated  to  be  +10  percent.  The 
theoretically  predicted  values  of  thrust  generally  fell  within  this 
+  10-percent  error  band.  The  theoretically  predicted  values  of  torque 
were  generally  slightly  higher  than  the  10-percent  error  band. 
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Since  the  experimental  data  were  not  considered  to  be 
sufficiently  accurate  and,  furthermore,  included  interference 
effects  (propeller -propeller  overlap,  fuselage  and  wing -tip  end 
effects)  not  included  in  the  present  mathematical  model,  no  definitive 
evaluation  of  the  theoretical  method  was  made.  Instead,  it  was 
recommended  that  accurate  performance  data  be  collected  on  an 
isolated  propeller  operating  in  hovering,  transitional,  and  axial 
flight  for  the  purpose  of  correlation  with  the  newly  developed  theory 
and  any  subsequent  theories. 

As  an  outgrowth  of  that  recommendation,  the  effort  reported 
herein  was  undertaken.  The  purpose  of  these  tests  was  to  obtain 
thrust,  torque,  and  inplane  forces  and  moments  on  an  isolated 
propeller  operating  in  hovering,  transitional,  and  axial  flight. 

Blade  strain  data  were  also  to  be  obtained.  The  parameters  which 
were  varied  were  free-stream  velocity  (  ),  propeller  speed  (CL  ), 

blade  angle  (fiK  ),  and  angle  which  the  thrust  axis  makeB  with  the 
free-stream  velocity  (r  ).  Correlation  with  theoretically  predicted 
propeller  performance  and  blade  stresses  was  to  be  performed. 


2 


DESCRIPTION  OF  THE  TEST  APPARATUS 


SMALL  ROTOR  TEST  FACILITY 

The  tests  were  conducted  on  a  test  bed  developed  by  CAL  for 
the  U.  S.  Army  in  1959.  The  facility, shown  in  Figure  1, consists  of  a 
modified  tractor -trailer  combination  capable  of  speeds  up  to  60 
miles  per  hour.  The  tractor  is  basically  a  standard  truck  chassis 
whose  body  was  modified  to  be  aerodynamically  streamlined.  At  the 
rear  of  the  tractor  is  a  rearward-facing  seat  for  the  test  observer 
and  the  remote  controls  for  the  test  conditions  and  data  recording 
on  the  trailer.  The  trailer  consists  of  a  large  flat  bed  from  which 
protrudes  a  tower  whose  height  is  adjustable  in  fixed  increments. 

Slung  beneath  the  trailer  bed  are:  two  motor -generator  sets  which 
supply  up  to  5,000  watts  of  power;  a  375-horsepower  automotive  engine 
which  supplies  power  to  drive  a  test  model  via  a  right-angle  gear 
drive  and  power  shafting  attached  to  the  tower;  water  and  gasoline 
supplies  for  the  engine;  a  radiator  and  an  oil  pump  for  cooling  of  the 
right-angle  gear  drive;  control  panels  for  operating  various  components 
of  the  facilities;  and  data  recording  equipment.  Mounted  at  the  top  of 
the  tower  is  a  pully  box  which  transfers  power  from  the  main  drive  shaft 
coming  up  from  the  engine  to  a  second  shaft  whose  centerline  is  coin¬ 
cident  with  the  centerline  of  the  tower.  The  entire  trailer  assembly 
is  mounted  on  an  air  suspension  system  designed  to  cushion  the  trailer 
from  road  roughness  and  to  maintain  the  trailer  bed  horizontal  at  all 
times.  For  a  more  detailed  description  of  the  basic  facility,  see 
Reference  3. 

MODIFICATIONS  TO  THE  SMALL  ROTOR  TEST  FACILITY 

In  order  to  allow  testing  of  an  isolated  propeller  on  the  above 
facility,  a  large  pod  was  designed  and  fabricated  (see  Figure  2).  The 
pod,  which  is  mounted  atop  the  tower,  houses  a  right-angle  gear  drive 
(similar  to  the  one  slung  beneath  the  trailer  bed)  which  transfers  rotational 
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power  from  the  vertical  shaft  of  the  pully  box  to  a  horizontal  shaft. 
This  shaft  is  supported  by  two  angular  contact  bearings  and  itself 
supports  a  slip-ring  assembly,  force  balance  unit, and  the  propeller. 
The  pod  has  a  counterbalance  attached  to  its  rear  face  to  insure  that 
the  center  of  gravity  of  the  entire  unit  is  located  on  the  centerline 
of  the  tower.  The  pod  assembly  is  mounted  on  a  turntable  which 
allows  rotation  in  a  horizontal  plane  through  360  degrees  in  5-degree 
increments. 

A  radiator  and  an  oil  pump  for  cooling  the  upper  right-angle 
drive  were  also  added, 

! 

I 

The  slip-ring  assembly  mounted  on  the  shaft  in  the  pod 
has  44  rings.  The  unit  is  supported  on  2  ball-bearing  rings, 
one  at  each  end.  The  silver  graphite,  chemopolarized  brushes  are 
mounted  on  readily  removable  brush  blocks  and  ride  on  coin-silver 
rings. 


The  propeller  set  which  was  tested  consisted  of  three  right- 
handed  blades.  These  blades  are  similar  to  the  ones  employed  on  the 
0.  6  scale  model  of  the  XC-142.  The  blade  planform  characteristics  are 
presented  in  Figure  3.  The  blades  were  mounted  in  a  standard  hub 
assembly  with  a  pitch-change  mechanism  powered  by  a  28-volt  dc 
supply.  Although  the  capability  for  changing  blade  pitch  angle  remotely 
could  have  been  incorporated,  it  was  felt  that  it  was  not  worth  the 
additional  effort  required.  Thus,  provisions  were  made  which  per¬ 
mitted  a  blade-angle  change  only  when  the  prop  was  stationary. 

MEASURING  DEVICES 
Rotating  Balance 

To  measure  the  instantaneous  values  of  the  shaft  loads,  a 
four -component  balance  system  designed  and  built  by  NASA  was  em¬ 
ployed.  This  unit  was  mounted  on  the  end  of  the  power  shaft,  and  the 
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prop  was  then  hung  from  the  opposite  end  of  the  balance.  The  balance 
rotated  with  the  prop  and  shaft.  Thus,  all  loads  from  the  prop  were 
transferred  to  the  shaft  via  the  balance. 


The  balance  measures  thrust  and  torque  directly.  However, 
the  character  of  the  balance  is  such  that  it  measures  only  a  compo¬ 
nent  of  the  inplane  force  and  moment  (i.  e.  ,  the  force  and  moment 
acting  in  the  plane  of  the  prop).  In  order  to  determine  the  magnitudes 
and  directions  of  the  total  inplane  force  and  moment  from  the  balance 
outputs,  these  loads  must  be  periodic.  Further  discussion  of  the 
balance  is  presented  in  subsequent  sections  and  in  Appendix  I. 

Strain  Rosettes 

The  instantaneous  strain  on  the  camber  side  of  the  blade  was 
measured  at  four  spanwise  stations  (see  Figure  4).  At  each  station, 
a  45-degree  strain  rosette  (three  strain  gages  mounted  at  45-degree 
angles  with  respect  to  one  another)  was  bonded  to  the  blade  surface. 

The  gages  were  SR-4,  type  FABR-50-35/S-1  3.  The  resistance  of  each 
gage  in  the  rosette  was  350  +  2.  5  ohms  and  each  had  a  gage  factor  of 
2.  1  3  +_  1  percent.  Each  gage  of  each  strain  rosette  was  wired  into  one  arm 
of  a  Wheatstone  bridge.  The  remaining  three  arms  of  the  bridge  were 
made  up  of  0.  5watt3,  350  ohms  +  0,  1-percent  resistors. 

Miscellaneous 

To  measure  wind  speed  and  direction,  two  wind-speed  indi¬ 
cators  and  one  wind-direction  indicator  were  mounted  at  the  front 
of  the  trailer  deck.  One  of  the  speed  indicators  was  connected  to  a 
tachometer  in  the  tractor  and  was  used  by  the  driver  to  set  the  desired 
test  forward  velocity.  The  second  wind-speed  indicator  signal  and  the 
wind-direction  indicator  signal  were  recorded  during  each  test  run  on  the 
oscillograph  records  for  data  reduction  purposes. 


5 


The  minimum  blade  pitch  angle  was  determined  by  the  low- 
pitch  stop  (a  mechanical  limit)  which  was  an  integral  part  of  the  hub 
assembly.  There  was  a  similar  high -pitch  stop.  Both  were  adjustable 
within  limits.  The  low-pitch  stop  was  set  at  its  minimum  value  which 
for  the  blades  tested  was  5  degrees  at  the  64.  4-percent  span.  The  high- 
pitch  stop  was  set  at  20  degrees  at  the  64.  4-percent  span.  The  inter¬ 
vening  span  was  divided  into  1 -degree  increments  from  5  degrees  to 
20  degrees.  Thus,  by  carefully  lining  up  a  pointed  indicator  mounted 
on  the  hub  with  scribe  marks  on  the  blade  shank,  the  blade  angle  at 
the  64.  4-percent  span  could  be  determined  to  within  about  +0.  5  degree. 


The  prop  speed  was  measured  using  an  electromagnetic  device 
mounted  in  the  pod  assembly.  The  completion  of  a  magnetic  circuit 
between  a  pin  and  a  plunger  once  per  revolution  generated  an  electrical 
pulse  which  was  recorded. 

RECORDING  EQUIPMENT 

All  recording  equipment  was  mounted  on  the  test  facility  (see 
Figure  1). 

Recording  was  performed  on  a  3-kilocycle  carrier  system  con¬ 
sisting  of  16  amplifiers,  2  oscillographs,  and  2  magazines.  Each  balance 
and  strain-gage  signal  was  fed  into  an  amplifier  and  then  into  the  oscillo¬ 
graph.  Signals  recorded  directly  on  the  oscillograph  were:  prop  speed, 
wind  velocity  and  direction,  a  100-cps  timing  trace,  and  a  synchronizing 
trace  which  allowed  correlation  in  time  between  the  two  oscillographs. 

The  recording  equipment  could  be  activated  remotely  from  the 
observer's  post  in  the  tractor. 

In  addition  to  the  above  recording,  certain  of  the  data  signals 
were  displayed  on  a  panel  in  the  observer's  post  for  the  purpose  of  con¬ 
trol  and  as  a  safety  check.  Figure  5  shows  this  panel.  Displayed  are 


engine  and  prop  speeds  for  the  purpose  of  control.  For  the  purpose  of 
safety,  the  upper  and  lower  gearbox  oil  temperature  as  well  as  prop 
thrust  and  torque  were  monitored  —  the  latter  so  that  maximum  safe 
prop  loads  would  never  be  exceeded,  and  the  former  to  warn  of  danger¬ 
ously  high  temperatures  in  the  right-angle  gear  drives. 

Also  located  in  the  observer's  post  was  the  control  which 
governed  propeller  speed. 


CALIBRATION  OF  THE  MEASURING  DEVICES 


CALI  BRA  TION  OF  THE  ROTATING  BALANCE 
Pretest  Calibration 

Special  attention  was  given  to  the  unusual  characteristics  of 
the  rotating  balance.  For  example,  care  was  taken  so  that  the  inplane 
force  and  moment  calibrations  were  made  along  the  most  sensitive 
balance  axes.  The  calibration  matrix  obtained  from  static  calibration 
of  the  balance  is  presented  in  Table  I.  The  coefficients  are  given  in 
terms  of  load  per  meter  unit.  Thrust  and  NF  coefficients  are  in  pounds 
per  meter  unit  .and  torque  and  PH  are  in  foot-pounds  per  meter  unit. 

The  meter  unit  is  defined  as  the  deflection  which  would  be  obtained  on 
a  measuring  device  (meter)  when  a  prespecified  resistance  is  placed 
in  parallel  across  a  prespecified  arm  of  one  of  the  balance  bridges. 

The  numbers  in  parentheses  in  Table  I  are  those  values  obtained  when 
the  balance  was  loaded  negatively.  and  pm  are  the  measured 

components  of  the  total  inplane  force  and  moment.)  For  a  definition  of 
positive  loads,  see  Figure  6. 


TABLE  I 

ROTATING  BALANCE  CALIBRATION  MATRIX 


THRUST 

TORQUE 

NF 

PM 

THRUST 

1.666536 

(1.6868710) 

0.031305 

(-0.000199) 

-0.02627 

(-0.015816) 

0.001836 

(-0.031632) 

TORQUE 

0.0000270 

(-0.000290) 

0.85573 

(0.857619) 

-0.01  5686 
(-0.007320) 

0.001435 

(0.004193) 

HF 

-0.0008430 

(-0.0008370) 

0.012272 

(0.013874) 

0.589132 

(0.592773) 

-0.060163 

(-0.0615910) 

PH 

-0.0082910 

(-0.008208) 

0.028247 

(0.029394) 

-0.255183 

(-0.249581) 

0.809439 

(0.800671) 
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The  value  of  the  resistance  required  to  produce  1,000 
meter  units  of  deflection  in  each  performance  channel  is  presented 
in  Table  II. 

Table  IT 

RESISTANCE  REQUIRED  TO  PROVIDE  1,000  METER  UNITS  OF  DEFLECTION 


CHANNEL 

RESISTANCE 

THRUST 

TORQUE 

NF 

PM 

1 .60  x  105  ohm* 
0.8>1  x  105  Ohm* 
1.22  x  105  ohms 
2.32  x  IQ*  ohm* 

When  the  above  resistances  were  placed  across  their  re¬ 
spective  bridge  arms  and  recorded  on  the  oscillograph  equipment 
described  previously,  the  resultant  trace  deflections  obtained  were 
(for  the  noted  attenuator  settings  on  their  amplifiers)  those  indicated 
in  Table  III. 


Table  III 

OSCILLOGRAPH  TRACE  DEFLECTIONS  FOR  PRESCRIBED  RESISTANCE 


CHANNEL 

ATTENUATION 

SETTING 

OSCILLOGRAPH 

TRACE  DEFLECTION 

THRUST 

100 

2.40  inches 

TORQUE 

200 

2,18  inches 

NF 

150 

2. 1 4  inches 

PM 

too 

1.78  Inches 
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The  attenuator  setting  represents  the  relative  increase  or 
decrease  of  the  output  signal  accomplished  by  the  amplifier  for  a 
given  input  signal.  This  could  be  set  at  prespecified  settings  from 
1  to  1, 000.  The  largest  increase  in  output  signal  was  obtained  at  a 
setting  of  "1". 

Thus,  to  compute  the  calibration  matrix  elements  in  terms 
of  load  per  inch  of  oscillograph  deflection,  consider,  for  example, 
the  thrust  due  to  thrust  element.  To  convert  1.  686536  pounds  of 
thrust  per  meter  unit  into  pounds  of  thrust  per  inch  of  oscillograph 
deflection  at  any  given  attenuation,  proceed  as  follows: 

aTr  '  (o,tt)(R7t)(0tt)(Gtt)  (1) 


where 

a.Tr  =  main  thrust  element  in  pounds  of  thrust  per  inch 
of  oscillograph  deflection 

a-rr  =  main  thrust  element  in  pounds  of  thrust  per  meter 
unit  —  1.  686536  pounds  per  meter  unit 

R-rr  =  number  of  meter  units  caused  by  a  prespecified 
resistance  in  parallel  with  a  prespecified  arm 

of  the  balance  bridge  — 1,000  meter  units 

- 5 - 

1 .  6  x  10  ohms 

Ott  -  ratio  of  prescribed  resistance  to  resultant  oscillo¬ 
graph  trace  deflection  at  a  given  amplifier  atten- 

5 

uation  —  1.6x10  ohms/2.  4  inches 

6-rr  =  ratio  of  new  attenuation  setting  to  attenuation  setting 
when  determining  Ott  — new  attenuation/ 1 00 
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Posttest  Calibration 


After  the  data  had  been  collected,  it  was  decided  to  perform 
a  calibration  of  the  balance  while  it  was  mounted  on  the  test  facility. 

The  purpose  was  simply  to  obtain  a  check  of  the  main  diagonal  elements 
of  the  matrix  of  Table  II.  The  calibration  of  the  torque,  NF,  and  PM 
components  was  done  only  statically^while  the  thrust  was  calibrated 

both  statically  and  with  the  propeller  shaft  rotating  at  approximately 
1200  rpm. 

The  static  calibration  of  torque,  NP ,  and  PM  yielded  the  same 
main  diagonal  coefficients  as  obtained  in  the  pretest  calibration  to 
within  5  percent  on  NF  and  PM  and  2  percent  on  torque. 

The  static  and  dynamic  calibrations  of  thrust  indicated  that 
a  large  hysteresis  loop  exists  in  the  thrust  load  versus  balance  output 
curve.  A  nondimensionalized  form  of  his  curve  is  presented  in 
Figure  7.  The  slope  obtained  while  loading  the  balance  is  approxi¬ 
mately  10  percent  higher  than  the  slope  obtained  upon  unloading  the 
balance.  This  hysteresis-loop  measurement  was  repeated  several 
times  with  the  balance  stationary  as  well  as  rotating.  The  mean  of 
the  slopes  obtained  in  these  posttest  calculations  agreed  to  within 
2  percent  with  the  slopes  obtained  in  the  pretest  static  calibration 
where  the  data  reduction  techniques  had  masked  the  hysteresis  loop. 

During  the  dynamic  calibration,  certain  other  anomalous 
behavior  of  the  balance  was  noted.  Specifically,  it  was  noted  that 

(i)  When  a  constant  inplane  force  and  moment 

existed,  the  response  in  the  thrust  component 
output  of  the  balance  due  to  these  inplane 
loads  varied  nonlinearly  with  the  azimuthal 
orientation  of  the  load  relative  to  the  blades. 

(ii)  A  similar,  but  smaller,  effect  was  observed  in 
the  indicated  torque  response  to  inplane  forces 
and  moments. 
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(iii)  The  response  in  the  NF  component  output  of 

this  balance,  due  to  a  constant  applied  inplane 
moment, also  varied  nonlinearly  with  azimuth 
as  did  the  response  in  the  PM  component  due  to 
a  constant  inplane  force. 

To  substantiate  these  observations,  a  calibration  of  the  thrust 
component  was  devised  in  which  a  fixed  normal  force  load  of  581  pounds 
was  hung  in  the  plane  of  the  propeller.  With  the  balance  in  a  known 
angular  position,  specified  as  zero,  a  pure  thrust  load  of  known  magni¬ 
tude  was  applied.  A  record  of  all  balance  components  was  made.  With¬ 
out  changing  any  of  the  loads,  the  balance  was  indexed  by  hand  through 
prespecified  angular  positions  and  records  were  taken  at  each  position. 
The  balance  was  returned  to  the  zero  angular  position  ,  the  next  thrust 
load  increment  was  added, and  the  procedure  was  repeated.  The  slope 
of  the  thrust  load  versus  the  thrust  deflection  was  a  constant  for  all 
azimuthal  positions  and  agreed  with  the  results  previously  obtained  to 
within  +2  percent.  However,  the  variation  of  the  thrust  response  due  to 
the  inplane  load  was  a  nonlinear  function  of  azimuthal  orientation  of  the 
balance.  This  is  shown  in  Figure  8  where  plots  of  T ,  NF,  and  pm 
responses  in  the  form  of  nondimensionalized  oscillograph  deflection 
versus  azimuthal  orientation  of  the  balance  are  presented.  Data  were 
recorded  every  20  degrees  of  azimuth  for  each  applied  thrust  load. 

The  deflection  in  each  balance  component  at  for  each  thrust  load 

was  subtracted  from  the  respective  deflections  at  each  of  the  other 
azimuthal  positions.  This  effectively  removed  the  deflections  due  to 
the  thrust  load.  The  trace  deflections  which  remained  were  due  to 
coupling  from  the  constant  inplane  force  and  moment.  Further,  since 
the  inplane  force  and  moments  were  constant,  the  coupling-induced 
deflections  all  collapsed  to  the  single  curves  shown  in  Figure  8  for 
all  thrust  loadings.  The  variation  in  torque  due  to  the  inplane  load 
was  negligibly  small.  Note  that  all  of  the  responses  are  sinusoidal. 
However,  the  Pm  component  is  approximately  180  degrees  out  of  phase 
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with  the  NF  component, and  the  r  component  is  approximately  130 
degrees  out  ot  phase  with  NF . 

From  the  above  calibrations  and  several  other  selective 
tests,  it  was  concluded  that 

(1  )  The  main  diagonal  terms  of  thrust  and  torque 
in  the  calibration  matrix  were  independent  of 
the  magnitude  of  the  applied  load  and  of  the 
azimuthal  orientation  of  the  balance. 

(ii)  For  any  given  azimuthal  position  of  the  balance, 
the  slope  of  the  applied  inplane  load  versus 
balance  output  for  that  channel  was  independent 
of  the  load.  Further,  the  angular  position  of 
maximum  balance  sensitivity  to  inplane  loads 
was  known,  and  the  corresponding  calibration 
constants  were  those  employed  in  the  calibration 
matrix  of  Table  I. 

(iii)  The  effects  of  inplane  loads  on  the  outputs  of  the 
thrust,  NF,  and  .^components  of  the  balance 
were  nonlinear  functions  of  the  azimuthal  position 
of  the  balance.  The  nonlinear  coupling  effect 
was  more  pronounced  between  hf  and  Fm  than 
between  either  Wf  and  T  or  PM  and  T  . 

To  define  the  magnitudes  of  the  nonlinear  coupling  terms 
accurately  would  have  required  considerable  additional  effort.  It  was 
decided  that  this  effort  was  unwarranted  at  this  time  in  view  of  the  +5 
percent  uncertainty  in  the  thrust  calibration  constant  due  to  the  hysteresis 
loop  and  the  belief  that,  even  if  the  coupling  terms  listed  in  Table  I 
were  in  error  by  100  percent,  their  effect  upon  the  evaluation  of  the 
mean  thrust  load  would  be  small.  Hence,  the  matrix  of  Table  I  was 
employed  to  reduce  the  experimental  data. 
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While  the  mean  loads  will  probably  not  be  much  affected  by 
the  errors  due  to  the  uncertain  magnitude  and  character  of  the  coupling 
terms,  the  fundamental  and  higher  harmonics  will  be  seriously  affected, 

CALIBRATION  OF  THE  STRAIN  GAGES 

A  test  beam  consisting  of  a  cantilevered  aluminum  bar  6.  25 
inches  by  1. 75  inches  by  0.  0625  inch  was  instrumented  with  a  strain 
rosette  of  the  same  type  mounted  on  the  prop.  The  test  rosette  was  wired 
into  a  bridge  network  similar  to  the  ones  employed  on  the  prop.  The 
output  from  the  bridge  was  connected  to  an  open  set  of  slip  rings  to 
allow  recording  of  the  signal  on  the  equipment  mounted  on  the  test 
facility.  The  beam  was  then  loaded,  incrementally,  with  a  tip-applied 
load, and  records  were  taken  at  each  loading.  The  test  was  repeated 
with  the  bridge  output  connected  to  a  strain  analyzer.  The  strain  of 
each  gage  at  each  applied  load  was  read  and  recorded. 

Plots  of  applied  load  versus  oscillograph  deflection  and  strain 
were  made.  From  these  plots,  a  plot  of  strain  versus  oscillograph 
deflection  was  obtained.  From  this  plot  was  obtained  the  strain  gage 
calibration  constant.  To  account  for  anticipated  differences  in  the 
electrical  characteristics  of  the  rosette  systems  mounted  on  the  blade, 
a  resistance  was  computed, which  when  placed  in  parallel  with  the  bridge 
arm  containing  a  strain  gage, resulted  in  a  1-inch  oscillograph  trace 
deflection.  The  shunt  resistance  was  then  placed  across  each  gage  of 
the  test  beam  and  the  blade.  The  test  beam  rosette  output  yielded  a 
1 -inch  trace  deflection  for  each  of  its  gages  with  an  error  of  approxi¬ 
mately  +0.  5  percent.  The  trace  deflection  obtained  from  each  prop 
strain  gage  thus  gave  a  measure  of  the  resistance  differences  in  each 
circuit  due  to  line,  slip  rings.and  amplifiers  and  could  be  used  to 
modify  the  calibration  constant  obtained  from  the  test  beam  strain 
rosette.  With  the  amplifier  attenuations  for  all  strain  gages  set  at 
"7",  the  resultant  calibration  constants  are  those  presented  in 
Table  IV. 
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Table  H 

STRAIN-GAGE  CALIBRATION  CONSTANTS 


STATION 

CALIBRATION  CONSTANT 

INCHES/INCH  OF  TRACE  DEFLECTION 

1  -  GAGE  1 

31.2485  *  10-6 

GAGE  2 

20.9065  *  10”6 

GAGE  3 

30.6589  x  10-6 

2  -  GAGE  1 

30.0911  x  10-6 

GAGE  2 

31.3154  x  10-# 

GAGE  3 

26.9324  x  10"6 

3  -  GAGE  1 

25.3016  x  10_# 

GAGE  2 

30.3724  x  I0-6 

GAGE  3 

62.3638  x  10"6 

4  -  GAGE  1 

29.3366  x  10"6 

GAGE  2 

28.5352  x  10-6 

GAGE  3 

32.0006  x  10"# 

------  ----- 
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CALIBRATION  OF  THE  BLADE  PITCH  ANGLE 

The  hub  with  the  blades  mounted  was  placed  on  a  level  refer¬ 
ence  plate.  A  reference  mark  was  made  on  each  blade  at  the  64.  4- 
percent  span.  Using  an  inclinometer,  the  blade  angle  at  the  low-pitch 
stop  position  was  measured.  It  was  found  to  be  5.  200  degrees, 

5.  317  degrees, and  5.  167  degrees  for  the  three  blades  measured. 

At  the  high-pitch  stop  position,  the  blade  angles  measured  were 
20.  083  degrees,  19.  700  degrees.and  20.  384  degrees.  To  determine 
how  well  intermediate  blade  angles  could  be  set,  the  blade  was  set 
using  the  pitch-change  mechanism  at  10  degrees  and  15  degrees  many 
times.  It  was  found  that  the  intermediate  angles  could  be  set  to  within 
+0.  5  degree  of  the  desired  values. 

CALIBRATION  OF  THE  WIND-SPEED  AND  -DIRECTION  INDICATORS 

The  wind-speed  and  -direction  indicators  which  were  recorded 
on  the  oscillograph  were  calibrated  to  check  the  expressions  supplied 
by  the  manufacturer . 

The  wind-direction  calibration  slope  used  in  data  reduction 
is  59.  0  degrees  per  inch  of  oscillograph  deflection.  The  associated 
error  is  +6  percent. 

The  wind-speed  calibration  supplied  by  the  manufacturer 
related  the  meter  output  in  cps  to  the  wind  speed  in  miles  per  hour, 

Vf  *  (3.  384381)  +  0.  627  (2) 

where 

Vf  -  the  indicated  velocity  in  mph 
u) -  the  measured  output  frequency  in  cps 

From  calibrations  performed  in  CAL’s  wind  tunnel,  Equa¬ 
tion  (2)  yields  velocities  approximately  2  percent  lower  than  actual, 
provided  the  wind  speed  is  greater  than  7  feet  per  second.  For 
wind  speeds  less  than  7  feet  per  second,  the  error  increases  rapidly. 
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TEST  PROCEDURES 


The  tests  were  conducted  on  the  main  northeast-southwest 
runway  at  the  Greater  Buffalo  International  Airport,  Buffalo,  New  York. 
Days  on  which  tests  could  be  conducted  were  limited  by  weather  and 
wind  conditions.  Tests  were  not  conducted  unless  the  local  wind  speed 
was  less  than  3  mph.  This  requirement  generally  forced  testing  to  be 
conducted  during  the  very  early  morning  hours. 

Recording  instrumentation  was  allowed  at  least  4  hours  to 
warm  up  prior  to  collecting  data.  This  was  required  to  insure  stable 
operation  of  the  amplifiers,  A  safety  checkout  procedure  was  executed 
prior  to  each  set  of  runs. 

Atmospheric  temperature  and  pressure  were  recorded  prior 
to  and  at  the  end  of  each  day's  running. 

Since  the  most  difficult  parameter  to  change  was  the  angle 
between  the  thrust  axis  and  the  free-stream  velocity  (r),  this  angle 
was  set  the  day  before  an  anticipated  test  date  and  generally  not  changed 
until  all  other  parameters  had  been  varied  over  their  entire  range. 

With  r  fixed,  the  blade  angle  (^)  was  set.  The  reference  blade  was 
placed  in  a  horizontal  position  and  all  amplifiers  were  balanced.  A 
zero  reference  record  was  taken.  The  propeller  was  then  brought  to 
a  preselected  speed  (/l  ).  Each  signal  trace  was  checked  to  insure 
that  it  was  neither  so  large  that  it  was  off  the  paper  nor  so  small  as 
to  make  its  analysis  difficult.  Any  necessary  changes  in  amplifier 
attenuations  were  made.  Preliminary  test  condition  data  were 
recorded;  i.  e.  ,  indicated  prop  speed,  attenuation  settings  of  all 
amplifiers,  atmospheric  pressure  and  temperature,  blade-angle 
setting,  etc.  With  the  prop  speed  at  /I  and  no  forward  speed,  a  data 
point  was  taken.  This  point  constituted  a  hover  point  (provided  the 
local  wind  velocities  were  negligibly  small).  Then,  holding  n  fixed, 
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the  test  facility  was  accelerated  to  a  predetermined  velocity  (  Vf  )  and 
maintained  at  that  speed  as  long  as  runway  length  permitted.  Records 
were  taken.  At  the  same  A  ,  data  at  several  other  Vf's  were  taken. 

Then,  with  T  and  fiR  held  at  the  same  values  and  the  test  facility 
stationary,  the  next  A  was  set  and  the  trace  levels  were  checked.  Pre¬ 
liminary  test  condition  data  were  recorded  as  was  the  hover  data  point. 
The  test  facility  was  again  accelerated  to  the  desired  Vf.  and  data 
were  taken. 


When  all  the  desired  prop  speed  variations  had  been  accom¬ 
plished,  fiR  was  changed  and  the  same  variations  in  A  and  Vf  were 
performed. 


The  last  variable  to  be  changed  was  r  ,  The  above  procedure 
was  repeated  for  fiR,  A,  and  Vf  . 


A  given  A  and/or  was  generally  repeatable  to  within 
+10  percent. 
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DATA  REDUCTION 


Each  performance  trace  was  identified  and  read  at  15 
equally  spaced  intervals  for  one  propeller  revolution.  The  values 
were  digitized  and  fed  into  a  digital  computer  program.  The  pro¬ 
gram  performed  the  operations  necessary  to  convert  the  trace 
deflections  to  loads,  to  remove  the  tare  values  due  to  prop  weight, 
and  to  print  out  the  results.  The  program  has  the  capability  of  har¬ 
monically  analyzing  the  performance  data.  However,  due  to  the  fact 
that  the  fundamental  and  higher  harmonics  of  the  performance  quan¬ 
tities  are  seriously  affected  by  the  uncertainty  associated  with  the 
magnitude  and  character  of  the  calibration  matrix  (Table  I),  the 
harmonically  analyzed  performance  data  were  not  presented.  A 
listing  of  the  program  is  presented  in  Appendix  II. 

The  strain  data  were  reduced  by  reading  the  mean  value  of 
each  strain  trace  and  then  feeding  this  digitized  value  to  a  digital 
computer  program.  The  program  converted  the  trace  deflections 
into  strains  experienced  by  each  gage  at  each  spanwise  station.  These 
strains  were  then  combined  to  yield  the  principal  strains  and  principal 
axes.  The  principal  stresses  were  then  computed.  The  program  is 
capable  of  harmonically  analyzing  the  strain  data  also.  The  strain 
data  reduction  program  is  included  in  the  listing  of  Appendix  II. 
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EXPERIMENTAL  RESULTS 


A  total  of  60  data  points  was  collected.  A  data  point  was  any 
combination  of  Vf  ,  fL  ,  t  ,  and  at  which  the  prop  performance 
and  blade  stresses  were  recorded.  The  ranges  of  parameters  covered 
are  presented  in  Table  V.  For  each  D.  listed,  data  were  obtained  at 
four  values  of  Vf  .  The  relatively  low  values  of  j&R  selected  were 
required  by  the  limitation  on  the  available  torque  from  the  test  facility 
power  plant.  The  maximum  torque  theoretically  available  was  470 
foot-pounds.  The  maximum  torque  actually  obtained  during  these 
tests  was  360  foot-pounds. 

The  forward  velocity  limits  were  determined  by  the  runway 
length  available  and  the  acceleration  limitation  of  the  test  facility. 
Thus,  the  maximum  speed  attained  was  approximately  50  miles 
per  hour. 

The  accuracy  to  which  the  experimental  performance  data  is 
known  varies  from  component  to  component.  Consider  first  the  thrust 
component.  It  is  believed  that  the  margin  of  error  associated  with 
the  thrust  measurement  is  approximately  +8  percent.  This  is  due 
to  (i)  the  +5  percent  uncertainty  in  the  thrust  calibration  constant 
associated  with  the  experimental  hysteresis  loop,  (ii)  the  questionable 
magnitude  of  the  coupling  terms  due  to  the  inplane  force  and  moment 
in  the  thrust  equation, and  finally,  (iii)  the  errors  associated  with 
reading  the  traces.  The  torque  is  known  to  +4  percent.  The  inplane 
force  and  moment  measurement  errors  are  estimated  to  lie  between 
10  percent  and  20  percent.  This  is  due  largely  to  the  fact  that,  gen¬ 
erally,  the  inplane  loads  and  moments  measured  were  small  compared 
to  the  correction  necessary  to  account  for  the  propeller  weight  com¬ 
ponent  and  the  large  uncertainty  associated  with  the  magnitude  of 
the  coupling  terms.  The  strain  measurements  are  believed  to  be 
accurate  to  +6  percent. 
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Table  1 

TEST  PARAMETER  VALUES 


t.  r  =  0° 

A  =  5° 

j 

CL  c  800  rpm 

1380  rpm 

1 800  rpm 

Vj  =;  0  mph 

13  mph 

27  mph 

45  mph 

2.  r  =  30° 

> 

II 

CJ1 

e 

/I  ~  800  rpm 

Vf  w  0  mph 

1380  rpm 

13  mph 

1800  rpm 

27  mph 

rpm 

45  mph 

A  8  '0° 

CL  si  600  rpm 

Vf^:  0  mph 

1240  rpm 

27  mph 

3.  ^  =  45° 

A  =5* 

/I  sr  800  rpm 

Vj ~  0  mph 

1380  rpm 

13  mph 

1800  rpm 

27  mph 

45  mph 

• 

o 

II 

< 

A  a;  800  rpm 

Vfs;  0  mph 

27  mph 

45  mph 

A  =  15* 

CL  «  800  rpm 

Vf  0  mph 

27  mph 

45  mph 

4.  r  =  65' 

A  =  5' 

Cl  ~  800  rpm 

Vf  ~  0  mph 

1360  rpm 

13  mph 

1800  rpm 

27  mph 

15  mph 

> 

II 

o 

• 

.ft  s;  1380  rpm 

Vf  =•  0  mph 

23  mph 
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In  addition  to  the  error  in  the  wind-speed  measurement 
discussed  previously,  there  existed  an  error  in  this  measurement 
due  to  flow  induced  by  the  propeller.  For  the  hovering  and  very  low 
forward  speeds,  the  velocity  induced  by  the  propeller  at  the  location 
of  the  anemometer  was  not  negligible.  A  correction  to  the  measured 
airspeed  based  on  the  induced  flow  field  predicted  in  Reference  4  was 
made.  At  speeds  in  excess  of  30  feet  per  second,  the  correction  was 
negligible  and  the  wind  speed  measured  was  used.  It  is  estimated  that 
for  speeds  greater  than  7  feet  per  second, the  error  in  the  measured 
wind  velocity  is  +7  percent.  At  speeds  less  than  7  feet  per  second, 
the  error  may  be  as  high  as  20  percent. 

Plotted  in  Figures  9  through  16  are  the  values  of  the  experi¬ 
mental  performance  characteristics  versus  free-stream  velocity  for 
several  values  of  r  and  fiK.  Indicated  on  each  curve  are  the  experi¬ 
mental  values,  the  associated  error  band, and  the  measured  propeller 
speed  (n.  ).  Curves  were  faired  through  experimental  points  for 
which  fl  did  not  deviate  by  more  than  +5  percent  from  the  mean. 

The  thrust  versus  curves  generally  indicate  that.as  the 
velocity  is  increased,  the  thrust  output  decreases.  This  appears  to 
be  true  for  any  fl  ,  fin  ,  or  r  with  the  exception  of  r  =  65  degrees 
(Figures  15  and  16).  For  this  r ,  the  curves  indicate  an  increase 
in  thrust  to  some  value  of  l£,and  thereafter, a  decrease  in  thrust  for 
further  increases  in  .  The  peak  or  maximum  thrust  depends  upon 
the  prop  speed.  Note  that  for  the  very  low  values  of  jfi.  and  relatively 
high  Vf  ,  the  thrust  curve  became  negative  for  all  values  of  r  except 
T  -  65  degrees.  For  approximately  the  same  A  ,  the  thrust  curve 
became  negative  at  higher  values  of  ^  as  r  increased.  The  thrust 
curves  moved  to  higher  values  almost  parallel  to  themselves  as  A 
was  increased. 

Increasing  forward  velocity  for  a  fixed  fl  did  not  seem  to 
affect  the  required  torque  (q)  over  the  velocity  range  tested.  The 
torque  curves  (Q  vs  Vj  )  moved  to  higher  values  almost  parallel  to 
themselves  as  ft  was  increased. 
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While  the  variation  of  both  the  7  and  Q  curves  was  smooth 
and  generally  decreased  slightly  as  Vf  increased  for  all  values  of 
n  ,  r  ,  and  j8r  ,  this  was  not  the  case  with  the  inplane  force  and 
moment  curves.  The  inplane  force  was  plotted  in  terms  of  its  com¬ 
ponents:  normal  force  (Mr)  and  side  force  ( SF  ).  The  inplane  moment 
was  plotted  in  terms  of  its  components:  pitching  moment  (PM)  and 
yawing  moment  (  YM ). 

In  Figure  9  (c-f),  the  variation  of  the  inplane  force  and  mo¬ 
ment  components  is  plotted  versus  for  jSR  =  5  degrees  and  T  -  0 
degrees.  This  corresponds  to  axial  flight.  For  the  configurations  tested, 
no  inplane  force  or  moment  was  expected  due  to  aerodynamic  loading. 
However,  variations  with  Vf  and  Cl  were  observed  in  all  components. 

The  components  observed  might  be  attributed  to: 

(i)  Aerodynamic  forces  and  moments  generated  by 
the  hub  and  pitch-change  mechanism. 

(ii)  Aerodynamic  interference  effects  from  the  pod 
and  tower. 

(iii)  Deviation  from  a  true  axial  flight  condition. 

The  characteristics  of  the  and  curves  continue  to  behave 
in  a  peculiar  fashion  as  V  is  increased  (Figures  10  through  16).  The 
PM  and  YM  curves,  however,  seem  to  indicate  definite  trends.  With 
the  exception  of  the  curves  for  Cl  1  380  rpm  (  r  =  30  degrees  and 
fin  -  5  degrees),  Cl  =»  1440  rpm,  and  Cl  -  880  rpm  (  r  =  45  degrees 
and  =  5  degrees),  the  PM  decreases  as  Vf  increases  for  all 
values  of  Cl  and  r.  The  same  comment  appears  to  be  applicable  for 
fin  -  10  degrees  and  PR  =  15  degrees  (see  Figures  11,  13,  14,  and  16). 
The  behavior  of  P/^with  IT  is  not  clear. 

With  the  exception  of  axial  flight,  the  behavior  of  Wwith  Vt 
for  all  17  ,  r,  and  is  to  increase  as  Vf  increases  to  some  value, 
and  then  decrease  as  Vf  is  further  increased.  For  any  given  r  other 
than  0  degree  (  pR  fixed  at  5  degrees)  and  Vf  ,  an  increase  in  Cl 
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increases  YM .  Further,  as  Z  is  increased,  the  YAY  is  increased  for 
any  given  Vf  and  n  ( fiR  again  fixed  at  5  degrees). 

From  the  data  of  Figures  9  through  16,  cross  plots  were  made 
to  indicate  trends  with  Ci  ,  7,  and  fi r  .  Figures  17  through  20  are  plots 
of  T  and  0,  versus  the  average  prop  speed  noted  on  Figures  9  through 
16,  The  plots  are  presented  for  jSe  -  5  degrees  only  and  for  various 
values  of  Vf  and  r  .  Again,  the  7  and  Q  are  generally  smooth,  well- 
behaved  curves  showing  an  increase  in  both  thrust  output  and  torque 
required  as  71  is  increased  for  any  specified  7  and  Vf  , 

Figure  21  is  a  plot  of  T  and  Q  versus  T  for  fix  -  5  degrees. 
Shown  is  the  variation  of  T  and  Q  for  two  prop  speeds  and  two  veloc¬ 
ities  representing  the  upper  and  lower  experimental  bounds  of  these 
variables.  It  is  noted  that  as  Z  is  increased  from  0  degree,  the 
thrust  increases  to  a  maximum  at  some  value  of  T  .  The  value  of  T 
at  which  this  maximum  occurs  is  dependent  upon  Zf  and  fl  .  As  7 
is  increased  further,  the  thrust  begins  to  decrease. 

The  variation  of  the  required  torque  is  more  complex.  As 
r  is  increased  to  approximately  15  degrees,  the  required  torque 
increases  to  a  maximum.  A  further  increase  in  r  to  approximately 
50  degrees  results  in  a  decrease  in  torque.  Finally,  indications  are 
that, as  T  is  increased  beyond  50  degrees,  the  torque  will  begin  to 
increase  again. 

For  a  7  =45  degrees,  Cl  -  822  rpm.and  Vf  =  40  feet  per 
second,  Figure  22  shows  that  both  T  and  Q  increase  at  an  increasing 
rate  as  ft#  is  increased  from  5  degrees  to  15  degrees. 

Because  of  the  low  blade  angles  at  which  these  tests  were  run 
and  because  negative  thrusts  were  obtained  for  several  test  conditions, 
it  was  conjectured  that  the  propeller  was  being  operated  in  a  range  in 
which  its  efficiency  was  very  low.  To  substantiate  this  conjecture,  the 
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results  from  the  axial  flight  case  were  plotted  as  CT  ,  CP ,  and/7  versus 
J  where 

CT  -  thrust  coefficient  =  T/pn  D  (3) 

CP  -  power  coefficient  -  2 rr  Q /p n1 D5  (4) 

y  -  efficiency  =  (CT/Cp)J  (5) 

J  -  advance  ratio  .s  \/^  /nD  (6) 

3 

p  -  density  —  slugs/ft 

D  -  prop  diameter  —  ft 

/?  =  prop  speed  —  cps 

Presented  in  Figure  23  are  CT  ,  Cp  ,  and  rj  versus  J  based  on 
the  data  collected.  From  the  rj  versus  J  curve,  it  is  apparent  that 
at Pf  -  5  degrees  the  prop  is  very  inefficient. 

Plotted  in  Figures  24  through  27  are  the  measured  longitudinal 
stresses  versus  span  for  the  various  Cl  's  and  Vps  at  which  tests  were 
conducted.  The  strain  measurements  from  which  these  stresses  were 
computed  were  measured  on  the  top  or  camber  side  of  the  blade.  Pos¬ 
itive  stress  indicates  tension.  As  expected,  it  was  found  that  the  stress 
level  usually  decreased  as  the  spanwise  position  increased.  The  stress 
level  at  a  given  spanwise  position  generally  increased  as  Cl  increased. 
Further,  for  a  given  il  as  Vf  increased,  the  stress  level  increased. 
This  is  attributed  to  the  fact  that  as  Vf  increased,  the  thrust,  or  lift 
per  blade,  decreased.  The  decreased  blade  load  reduces  the  com¬ 
pression  load  induced  in  the  fibers  on  the  blade  camber  (atop)  surface. 
Thus,  the  rotation-induced  tensile  load  is  reduced  by  a  smaller  amount 
and  the  entire  longitudinal  stress  increased. 
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Figure  24  presents  the  results  for  T  =  0  degree.  For 
Cl  as  850  rpm,  the  load  distribution  is  such  that  the  observed 
longitudinal  stress  becomes  negative  at  the  3.  5-foot  radius  for  all 
forward  velocities.  (The  two  curves  designated  by  7  =  0  degree, 

-  5  degrees,  A  =  873  rpm,  and  Vf  -  16  miles  per  hour  were 
run  under  the  same  test  conditions  to  determine  repeatability.  ) 

The  reason  for  this  peculiar  distribution  of  the  stress  curve  is 
unknown.  As  Cl  is  increased,  the  tensile  stress  curves  become 
smoother  and  are  positive  at  all  spanwise  locations  except  for  the 
case  at  Cl  =.  1500  rpm  and  Vf  -  46  miles  per  hour. 

A  similar  behavior  is  observed  in  Figure  25  where  the  pri¬ 
mary  change  has  been  in  r  (  r  is  now  30  degrees).  The  behavior  of 
the  stress  curve  with  JT  and  is  substantially  the  same  as  at  T  =0 
degree.  An  increase  in  fiR  from  5  degrees  to  10  degrees  results  in  a 
marked  decrease  in  stress  level  for  the  same  A  and  .  This  is 
again  attributed  to  the  increased  load  generating  increased  compression 
loads  due  to  bending  in  the  upper  fibers  of  the  blade.  Note  that  the 
unusual  behavior  at  Cl  =  850  rpm  disappeared  at  -  10  degrees. 

The  results  for  T  -  45  degrees  are  presented  in  Figure  26. 

For  -  5  degrees,  the  effect  of  V/  on  the  stress  was  not  as  pro¬ 
nounced  as  for  Z  =0  degree  and  30  degrees.  However,  as 
increased  to  10  degrees  and  15  degrees,  the  effect  of  became 
larger.  Note  again  the  marked  decrease  in  blade  stress  on  the  camber 
face  as  is  increased. 

Figure  27  presents  the  results  for  r  =  65  degrees.  The 
free-stream  velocity  has  almost  no  effect  on  the  stress  level  for 
=  5  degrees  at  this  r  . 
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COMPARISONS  OF  THE  EXPERIMENTAL  DATA  WITH 
THEORETICAL  PREDICTIONS  AND  OTHER  EXPERIMENTAL  DATA 


An  attempt  was  made  to  compare  the  data  obtained  in  the 
current  tests  with  the  results  presented  in  Reference  5.  The  tests 
in  Reference  5  were  performed  in  the  NASA  Ames  40-foot  by  80-foot 
wind  tunnel  on  three  propellers  having  three  blades  each.  The  plan- 
forms  were  different  from  one  another  as  well  as  being  different 
from  the  blade  planform  used  in  the  current  tests.  However,  since 
the  trends  reported  were  much  the  same  for  all  the  propellers  tested 
in  Reference  5,  it  was  felt  that  the  propeller  tests  of  this  report  should 
yield  similar  trends. 

The  results  of  Reference  5  were  presented  as  plots  of  Cr  and 
Cp  versus  /'  at  various  A*  and  r  .  The  parameter  j'was  defined  as 

/'«  Vf  cos  T/nD  (Modified  Advance  Ratio),  (7) 

It  was  not  noted  whether  the  j'  variation  was  obtained  by  varying 
Vf  and  holding  Ci  constant,  or  by  varying  A  and  holding  Vf  con¬ 
stant,  or  by  allowing  both  Vf  and  CL  to  vary.  Cross  plots  were 
made  of  Cr  and  CP  versus  r  at  constant  values  of  /' .  One  of 
the  conclusions  drawn  was  that  "  Cr  and  Cp  for  given  values  of 
/Bg  and  /'  are  nearly  constant  in  magnitude  over  a  large  range  of 
r  and  that  range  diminishes  with  increasing  /'". 

The  data  presented  in  Figures  9,  10,  12, and  15  were  put  in 
coefficient  (Cr,  Cp  )  form  and  /'was  computed  for  each  point.  The 
curves  are  presented  in  Figure  28.  The  solid  curves  were  obtained 
holding  A  constant;  the  dotted  curves  were  obtained  holding  Vf  con¬ 
stant.  Thus,  the  character  of  the  Cr  and  Cp  versus  /'curves  appar¬ 
ently  depends  upon  the  manner  in  which  the  data  were  generated. 
Further,  if  cross  plots  of  CT  and  Cp  versus  r  at  any  specified  J' 
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were  made  from  either  the  i /f  =  constant  or  ft  =  constant  curves, 
the  curves  would  not  be  constant  with  r  over  any  range  of  T  .  These 
results  seem  to  be  in  contradiction  to  the  results  obtained  in  Refer¬ 
ence  5.  No  reasons  for  these  results  could  be  found  in  the  data  collected. 

Thus,  from  the  data  collected  by  CAL,  it  would  appear 
that  s'  is  not  an  independent  parameter  in  the  problem  of  a  propeller 
operating  in  transitional  flight. 

The  theoretical  computer  program  developed  by  CAL  (Refer¬ 
ence  1)  was  run  for  five  different  cases  for  the  purpose  of  corre¬ 
lation.  The  variables  associated  with  said  cases  are  presented  in 
Table  VI.  The  mean  values  of  the  performance  for  each  case  were 
plotted  in  Figures  9,  11.  12,  14,  and  16.  The  theoretical  points 
are  indicated  by  a  symbol  with  a  flag. 


Table  IT 

VALUES  OF  PARAMETERS  USED  IN  THEORETICAL  CALCULATIONS 


IDENTIFICATION 

E9 

A 

deg 

SI 

rpm 

H 

CASE  1 

0 

5 

1875 

50 

CASE  2 

30 

. 

10 

1277 

32 

CASE  3 

45 

5 

819 

31 

CASE  4 

45 

15 

634 

29 

CASE  5 

65 

10 

1428 

37 

The  axial  flight  results  of  Figure  9  indicate  that  the  predicted 
T  is  111  percent  too  high  while  the  predicted  Q  is  27.  9  percent  too  low. 
The  reason  for  this  very  large  discrepancy  in  thrust  is  attributed  to 
the  fact  that  the  propeller  was  operating  very  nearly  in  a  windmill 
brake  state.  The  theoretical  results  indicated  that  the  outboard 
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20  percent  of  the  blade  was  operating  at  negative  values  of  lift. 

This  operating  condition  causes  computational  errors  in  the  theory. 
First,  the  lifts  which  result  from  negative  effective  angles  of  attack 
are  in  error.  In  the  negative  angle-of-attack  range,  the  two- 
dimensional  airfoil  characteristics  were  programmed  to  be  simply 
straight-line  extrapolations  from  their  values  at  zero  angle  of  attack. 
This  approximation  becomes  progressively  poorer  as  the  angle  of 
attack  becomes  more  and  more  negative.  The  error  is  such  as  to 
predict  a  larger  local  lift  so  long  as  the  lift  is  positive.  When  the 
lift  becomes  negative,  the  error  is  such  as  to  predict  a  smaller  local 
lift.  Second,  the  fact  that  the  distribution  of  the  bound  vorticity  is 
such  that  it  passes  through  zero  three  times  (once  at  tip,  root,  and 
one  intermediate  spanwise  location)  presents  difficulties  in  assigning 
the  proper  strength  and  position  to  the  trailing  tip  and  root  vortices 
employed  in  the  mathematical  model  of  the  rolled-up  wake.  Indeed, 
a  more  reasonable  rolled-up  wake  representation  for  a  blade  having 
such  a  distribution  would  be  two  pairs  of  vortices. 

The  inplane  forces  and  moments  predicted  by  the  theory  for 
the  axial  flight  case  are  zero.  As  noted  before,  the  cause  of  the 
experimentally  obtained  forces  and  moments  is  unknown. 

The  transition  Case  2  results  are  plotted  in  Figure  11.  The 
predicted  thrust  is  15  percent  higher  and  the  torque  is  13  percent 
higher  than  the  corresponding  experimental  values.  The  values  of 
NF ,  SF  ,  and  PM  are  underpredicted  while  YM  is  overpredicted. 

In  Figure  12,  the  results  for  transition  Case  3  are  presented. 
The  thrust  was  overpredicted  by  approximately  17  percent  and  the 
torque  underpredicted  by  25  percent.  For  this  case,  the  theoretical 
a/f  and  PM  are  less  than  the  corresponding  experimental  values,  while 
the  theoretical  SF  and  YM  are  greater  than  the  experimental  values. 


29 


For  transition  Case  5  (Figure  16),  the  theoretical  thrust  is 
20  percent  higher  and  the  theoretical  torque  is  approximately  40  percent 
higher.  Again,  correlation  with  the  inplane  moments  and  forces  is  poor. 

During  the  theoretical  investigation  of  Case  5,  it  was  dis¬ 
covered  that  the  magnitude  and  position  of  the  tip  vortex  representing 
the  rolled-up  wake  played  a  significant  role  in  the  prediction  of  the 
thrust  and  torque.  This  was  found  to  be  so  even  though  the  wake  model 
was  such  that  the  effect  of  a  mesh  of  trailing  and  shed  vorticity  was 
included  in  the  computation  for  180  degrees  of  blade  azimuth  prior  to 
rolling  up  into  a  concentrated  tip  and  root  vortex.  It  was  found  that 
increasing  the  magnitude  of  the  rolled-up  tip  vortex  by  20  percent  and 
shifting  its  spanwise  position  in  the  wake  from  the  93-percent  span  to 
the  50-percent  span  caused  a  10-percent  reduction  in  the  predicted 
thrust  and  a  5-percent  reduction  in  the  torque. 

In  these  transition  cases,  it  was  noted  that  the  predicted 
effective  angle  of  attack  ( ae  )  was  generally  very  small  and  resulted 
from  the  differences  between  the  built-in  twist,  j5  ,  the  induced  angle 
of  attack,  .  and  the  live  twist,  9  ;  i.  e.  ,  at  any  station 

ag  *  -  ctjr  *  j  •  station  number.  (8) 

For  those  cases  in  which  JS  ss  az  ,  the  magnitude  of  6  plays 
a  very  important  role  in  determining  the  blade  loading.  Thus,  a 
relatively  small  error  in  the  predicted  &  distribution  could  result 
in  a  relatively  large  error  in  the  predicted  propeller  performance. 

In  an  attempt  to  support  the  above  arguments,  calculations 
for  Case  4  were  performed.  This  case  was  selected  because  it  was 
believed  that  the  final  effective  angles  of  attack  would  be  sufficiently 
large,  positive  quantities  that  the  errors  associated  with  both  nega¬ 
tive  angles  and  prediction  of  the  torsional  deformation  would  be 
minimized.  The  results  are  presented  in  Figure  14.  Agreement 
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between  the  theoretical  and  experimental  values  of  most  of  the 
performance  quantities  appears  to  be  much  improved  over  the 
previously  discussed  cases.  The  thrust  predicted  is  within  the 
estimated  experimental  error  associated  with  the  measurement. 
The  torque  is  overpredicted  by  approximately  23  percent. 

The  poor  correlation  obtained  for  the  previously  discuss 
cases  may  also  be  attributable  to  the  following: 

(i)  Inaccuracies  in  setting  the  blade  angle,  fiK  . 
Calculations  indicate  that,  at  fin  -  5  degrees, 
a  10-percent  reduction  in  fin  results  in  a 
21 -percent  reduction  in  T  and  a  15-percent 
reduction  in  (J  .  It  is  believed  that  similar 
variations  would  be  obtained  at  higher  fix's , 

(ii)  The  use  of  two-dimensional  airfoil  data  in 
the  analysis.  No  experimental  proof  is  yet 
available  to  answer  the  fundamental  question 
of  whether  the  use  of  two-dimensional  airfoil 
section  data  is  reasonable  for  the  propeller. 

(iii)  The  possibility  that  the  errors  in  the  experimental 
data  presented  are  larger  than  estimated. 

(iv)  In  a  digital  program  of  the  complexity  of  the  one 
developed  in  Reference  1  and  employed  here  for 
the  purpose  of  correlation,  the  possibility  for 
program  errors  always  exists.  Though  the  pro¬ 
gram  has  been  thoroughly  checked,  the  unusual 
blade  load  distributions  predicted  for  the  cases 
computed  may  have  caused  the  program  to  operate 
in  an  erroneous  mode  heretofore  undiscovered. 


CONCLUSIONS 


Based  on  the  results  obtained  In  this  study,  the  following 
conclusions  were  drawn: 

(i)  For  the  range  of  parameters  tested,  the 
quantity  V,  co*  r/nO  is  not  an  independent 
parameter  of  a  propeller  operating  in  transitional 

flight. 

(ii)  As  r  is  increased  from  0  degree,  the  thrust 
increases  to  a  maximum  at  some  7  ,  dependent 
upon  ^  and  CL  ,  and  then  decreases  as  Z  is 
increased  further. 

(ui)  As  r  is  increased  from  0  degree,  the  required 
torque  increases  until  r  reaches  approximately 
15  degrees.  A  further  increase  in  r  to  approxi¬ 
mately  50  degrees  results  in  a  decrease  in  torque. 
The  required  torque  begins  to  increase  again  as 
r  is  increased  beyond  50  degrees. 

(iv)  The  effect  of  increasing  forward  velocity,  all 

other  parameters  held  constant,  is  to  decrease 
the  thrust  output  while  the  required  torque  remains 
essentially  constant. 

(v)  The  thrust  and  torque  both  increase  with  an 
increase  in  prop  speed. 

(vi)  Both  the  thrust  and  torque  increase  as  the  blade 
angle  of  attack  is  increaseu. 


(vii)  When  the  propeller  was  operating  in  a  flight 
condition  for  which  the  theory  was  developed, 
correlation  between  theory  and  experiment  was 
good. 

(viii)  Accurate  prediction  of  the  blade  torsional  deflection 
plays  a  very  important  role  in  determining  the  blade 
loadings  accurately. 

(ix)  The  position  and  magnitude  of  the  rolled-up  tip 

vortex  must  be  predicted  accurately  in  order  that 
the  thrust  and  torque  may  be  predicted  to  within 
the  accuracies  required  for  V/STOL-type  propellers 
throughout  the  entire  flight  range. 
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Figure  6.  SCHEMATIC  OF  PROPELLER  PLANE  AND  POSITIVE  SHAFT  LOADS 
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APPENDIX  I 


DERIVATION  OF  THE  INPLANE  FORCE  AND  MOMENT 
FROM  DATA  OBTAINED  ON  THE  ROTATING  BALANCE 


The  rotating  balance  designed  and  built  by  NASA  has  the 
characteristic  that  it  responds  only  to  a  component  of  inplane  force 
and  moment.  For  example,  consider  a  rectangular  beam  canti¬ 
levered  at  one  end  (see  sketch). 


The  beam  would  deflect  most  due  to  a  load  applied  along  the  x  axis. 
Hence,  the  x.  axis  would  be  referred  to  as  the  sensitive  axis.  In  the 
analogous  situation  in  the  balance,  this  axis  is  fixed  in  the  balance  and 
rotates  with  it.  The  balance  is  instrumented  to  respond  only  to  loads 
along  this  axis.  Thus,  at  any  instant  of  time,  the  balance  will  put  out 
a  signal  which  is  proportional  to  the  component  of  the  inplane  force 
(or  moment)  along  the  sensitive  axis.  Consider  Figure  29  where  X0  , 

Y0  ,  Z0  is  a  coordinate  system  fixed  in  the  prop  plane  such  that  the 
za  axis  lies  along  the  axis  of  rotation.  Then,  the  normal  force  com¬ 
ponent  (  a if)  and  the  yawing  moment  component  (  YM)  are  defined  to  be 
components  of  the  inplane  force  and  moment  in  the  X0  direction, and 
the  side  force  (■S F  )  and  pitching  moment  (PM)  are  defined  as  the  com¬ 
ponents  in  the  Y0  direction. 

The  angular  orientation  of  the  sensitive  axis  is  defined  by 
measured  from  the  X0  axis. 


71 


t  R  be  the  vector  representation  of  the  total  inplane 
force  (moment)  when  i?  •  a, 

Tj,Z  be  the  unit  vectors  along  Xe,Yot  and  "sensitive" 
axes  .respectively . 

oc  be  the  angle  between  the  sensitive  axis  and  the 
Xe  axis, 

be  the  angle  between  R  and  X0  . 


Figure  29.  SCHEMATIC  OF  RESULTANT  INPLANE  FORCE. 


The  component  of  R  along  the  sensitive  axis  is  given  as 


B  -  R  ■  ir. 


(9) 


But,  for  the  inplane  force,  say, 

r  =  nfT  +  sfT 

and 

v’  *  ( cos  oe. )  i  +  ( sin  oc)  j .. 


Thus, 


fl  -  NF  cos  a  +  5Fjin  <X  . 


(10) 


(ID 


(12) 
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For  a  three-bladed  isolated  propeller  in  steady  flight,  the 
loads  repeat  three  times  per  prop  revolution.  Hence,  every  120  degrees, 
ft  has  the  same  magnitude  and  orientation  as  shown  in  Figure  29, 
but  the  angular  position  of  the  sensitive  axis  has  shifted  by  120  degrees. 
The  equations  describing  this  condition  are 

B,  -  NF  cos  Oc  +  SF sir>  oc.  (13) 

Bz  -  -NF sin  (30  +  (X)  +  5F  cos  (30  +  oc)  (14) 


where 

Bf  =  the  balance  reading  at  <p  ‘  os 
3Z  =  the  balance  reading  at  V-  or+120- 

From  the  set  of  Equations  (13)  and  (14),  the  unknowns  NF  and 
■SF  may  be  solved  to  obtain 


A IF 


SF 


YJ 

2 

/T 


B,  cos (30  ■F  ex)  —  sin  Oe  J 
Bz  COS  Os  +  Bf  sin  (30  +  Oc)  j. 


These  values  of  A/^and  5F exist  every  120  degrees. 


(15) 

(16) 


For  the  balance  system  in  hand,  the  quantities  3i  and  Bz  repre¬ 
sent  the  output  signal  denoted  NF  at  two  azimuthal  positions,  oc  and 
a +120  degrees.  Thus,  by  knowing  the  angular  position  of  the  sensitive 
axis  with  respect  to  the  X0  axis,  the  NF  signal  may  be  sampled  in  such 
a  fashion  as  to  yield  the  azimuthal  history  of  NF  and  SF. 

The  same  comments  apply  for  the  pitching  and  yawing  moments. 
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APPENDIX  II 

LISTING  OF  THE  DATA  REDUCTION  PROGRAM 


INPUTS  REQUIRED 

Card  1  1CHK  0  -  Does  not  print  uncor¬ 

rected  performance  results 

1  -  Prints  uncorrected  results 


Card  2  EYM 

PR 


Young's  modulus  (p.  s.  i.  ) 
Poisson's  ratio 


PSIB  Angular  position  of  balance 

sensitive  axis  with  respect 
to  reference  blade  (degrees) 

DELPHI  Azimuthal  increment  at  which 

data  is  sampled  (degrees) 


Card  3  RM(I),  I  *  Radian  position  from  prop 

center  of  rotation  to  each 
strain-gage  location  (feet) 

Cards  4  to  7  PCC(l,J),  J*  f ,  V;  I-f,  +  Elements  of  performance 

calibration  matrix  (load 
per  inch  oscillograph 
deflection) 
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,  no  - » 


Card  8 


Cards  9  to  12  5CC  (1.  J),  1*1,3-, 
J-  1.4 


Cards  13  to  16  5ZD{J,K),  K  -  1.3 ; 

r ■  /.* 


Performance  Case 

Card  17  Case  10(1),  J*  1,8 

Type 

NPTS 

SCT  F 


YsAlN 

DY 


*nr»T./r"’r;  ’r’WHWB 


Distance  from  reference 
line  of  record  to  no-load 
trace  for  each  performance 
quantity  (inches) 


Elements  of  strain  cali¬ 
bration  matrix  (10  ^  inches 
per  inch  per  inch  oscillograph 
deflection) 


Distance  from  reference  line 
on  record  to  no-load  trace  for 
each  strain  quantity  (inches) 


Identification  of  cases  to  be 
analyzed 

P  -  Indicates  performance 
data  follows 

-  Indicates  strain  data 
follows 

Number  of  increments«at 
which  data  is  read 

F  -  Plot  unsealed  data 
T  -  Plot  scaled  data 

Minimum  Y  value  to  be  plotted 

Scale  of  Y  axis  per  inch 


■ 

si 


14 

% 
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♦Card  18 


PA  (I),  I  =  /.  y 


Attenuation  ratio  for 
performance 


Card  18 


Cards  19  to  22 


Cards  19  to  30 


*  In  performance, 


5  A  (It  J)  ,1  *  1,5 ;  Attenuation  ratio  for  strains 

/.  ¥■ 


PD(J,  J),  1-1,  MPTS;  Performance  deflections 
T  -  1 ,  V  from  reference  trace  (inches) 


50(1',  1,  <),  1*  1,  A IPPS)  Strain  deflections  from 
Km  1,3  ;  J*  I ,  V  reference  trace  (inches) 


first  case  should  be  gravity  tare  corrections. 
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. r> err?*-* ■• 


c 


* ihftc  ssvtol  l  ist.oeck.rff 

DIMENSION  seen, 4) , PCC ( 4, 4 1 ,C ASE I D I  8)  , PA (41  ,  POOL  ,4  I  ,T  I  30  I  .OH*  I  , 
1FMBI  30,2>,PSIR(  31>,NF  (3r),SF(3CI,PMnCI  ,YMI31|  ,SA(3,4|,S0(3n,4,  M  . 
2E(3)»TS(30,4),SS(30,4) 

DIMENSION  PZUI4), S/DI  4, 3) ,ASIGN(4,3) , TO (15| ,001 1 5  I , NFGI 151 , 

1 SFO( 15),YMG( 15) ,PMO(  15) ,RM|4) 

DIMFNSION  PFRF ( 41 ,  STRSS(4) 

RFAl  NF,NFG 
LOGICAL  SC  TF 

CAT  A  AS! ON/- 1.0,1. 0,1. r  . 1 .. , 1 . ( »- 1 .C , 1 . r  » -1 ..» ♦ l .., ).<  , l.C ,-1.0/ 
OATA  P/1HP/.NCASE/-/ 

DATA  PERF/24HT  0  NF  PM  / 

OATA  STR5S/24HR 1  R2  R3  R4  / 

CALL  OVOCHK(l) 

READ  (5,11  ICHK 

1  FOPMAT  (61 121 

READ  (5,21  EYM,PR,PSIH,DELPSI 

2  FOPMAT  (4E12.5I 

READ  (5,2)  IKM( I ) ,  I  *  1 , 4) 

READ  (5,21  ( (PCC( I , J I ,  J *  1 ,4 1 ,  t *1 , 4 ) 

READ  (5,21  (PZO( 11,1*1,4) 

READ  (5,50)  ( I SCC (I,J),I*),?),J*1,4) 

50  FORMAT  (3E12.5) 

READ  (5,51)  ((SZD(J,KI,K«1,3),J*1,4) 

WRITE  (6,31 

3  FORMAT  (5CH1EXPERIMENTAL  DATA  REDUCTION  PROGRAM  -  A. TRFNKA/SS/ ) 
WRITE  ( 6, 9G  I  ICHK 

90  FORMAT  I 6HL ICHK- , I  2 ) 

WRITE  (6,14)  EYM,PR,PSlB,DElPS! 

14  FORMAT  ( 5H0EYM-, E12.5,2X, 3HPR*, El 2.  5*2X,5HPSIB**E12.5*2X* 
17HDFLPS I «, El 2. 5  I 
WRITE  (6,55)  (RM(l),I*l,4) 

55  FORMAT  ( 3H0RM/1X,4E15.5) 

WRITE  (6,5)  ( (PCC (  I  »JI *  J*1 , 4 ) ,1*1,4) 

5  FORMAT  ( 4H0PCC/I 1X.4E15.5)) 

WRITE  (6,511  (PZO(  1 1,1*1,41 

51  FORMAT  (16HC.PZD  (  T,Q  ,NF  ,  PM)  /  l  X,4F  7. 2) 

WRITE  (6,4) 

4  FORMAT  (4H0SCCI 
DO  92  J ■  1, 4 

WRITE  (6,93)  STRSS(J), (SCCII, J), 1-1,3) 

93  FORMAT  ( IX, A2/1 X, 3E 15 . 5  ) 

92  CONTINUE 

WRITE  (6,52)  ( (SZOI J,K) ,K*1,3)  ,J*1 ,4) 

52  FORMAT  ( 18H0SZ0  (R 1 ,R2,R 3,R4) /IX, 4( 3F6. 2 , 3X) > 

YMF«EYM/(1.0-PR**2> 

RPD*0. 017453292 
PSIBR»PSIB*RPD 
OELPSR»OELPSl*RPO 
KPTS-36C .O/DELPS! 

DO  40  I  * 1,KPTS 

PS1R(  I) -FLOAT! l-ll*DELPSR 

TGI  I  1*0.0 

QG( II-O.O 
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NFG(  n-o.o 
SFG(  D-0.0 
PMG( I  1-0.0 
YMGI II-O.C 
40  CONTINUE 

10C0  READ  (5,6)  (CASE !0(  I  I  ,  I • 1 , S) ,  TYPE  »NPT  S*  SCTP »  YMIN.OY 

6  FORMAT  (8A6,A1,I3,L4,2F12.5I 

WRITE  (6,7)  (CASE I0«  I  I  •  I  *  1 • 81 ,  TYPE ,NPT S, SC  1 F , YMJN ,DY 

7  FORMAT  (lHI,  8A6/ IX , 5HT YPf « , A1 ,2X , 5HNPTS- t 3 , 2 X , SHSCTF- , L? ,  2X , 
15HYMIN»,E12.5,2X,3HDY«,E12.5> 

IF  (TYPE.EO.P)  GO  TO  IOC 
GO  TO  200 

IOC  RE AO  (5,21  4 PA( I), 1-1,4) 

WRITE  (6,8)  (PA( I ),  1-1,41 

8  FORMAT  ( 3HOPA/1X.4E15.5) 

WRITE  (6,9) 

9  FORMAT  ( 3HOPO I 

00  60  J-1,4 

READ  (5,27)  ( PD( I , J  I ,  I - 1 .NPTS) 

27  FORMAT  ( I5F4.2) 

WRITE  (6,61)  PERF(J),(PD(I,J),I-1,NPTS) 

61  FORMAT  (  lX,A2/( 1X.15F7.2)) 

6C  CONTINUE  . 

12  00  15  1-1, NPTS 

T(I)  •  0.0 
Q(I)  -  O.C 
DO  11  J-1,4 

PAPO  «  PA(J)*(P/D(JI-PO(I,J)) 

T(  I)  «  T(l)»PCC(l,JI*PAPO 
0(1)  ■  Q1I )*PCC(2,J)*PAP0 
11  CONTINUE 

00  16  L«l,2 

FMB( I, L 1-0.0 
00  17  J-1,4 

PAPO  -  PA(J)*(PZO( JI-P0(I,JI) 

FMB( (,L )-FMB( I,L)*PCC(l*2, J)*PAPD 
17  CONTINUE 
16  CONTINUE 
15  CONTINUE 

IF  (  ICHK.EQ.O)  GO  TO  80 
WRITE  (6,81) 

81  FORMAT  ( 24H1UNC0RRECTED  PERFORMANCE) 

WRITE  (6,25) 

60  00  20  I-l.NPTS 

ALPHA1-PS  IR (  D-PSIBR 
SINAl-SINI ALPHA  1 ) 

C0SA1*C0S( ALPHA 1 ) 

S IN A 2- ( 0.8660254* SINAI  .  5*C0SAl 1 
COS A 2- ( C . 8660254*COSA  1-0. 5*S INA1 1 
IT2-M0D1 I*4,NPTS)*1 
00  21  L-1,2 

6T1-FM6( I , L ) 

BT2-FMB ( I T  2 , L ) 

GO  TO  ( 23,24) ,L 

23  NF( I )»l .1547*(BT1*C0SA2-6T2*S  INAl ) 
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SF {  I  )*1,1547*(BT2*CQSA1*BT1*S INA2 I 
GO  TO  21 

24  PM<  I  l*1.1547*JBT2*CaSAl*BTl*SINA2» 

YM  (  f  )«  l.  1 547*  <BT>COSA2-BT2*  SINAI) 

IF  IICHK.EO.il  WkI’E  <6.26)  < I  ,  T( 1 1 ,Q (I) ,NF  (  I  I , SF < 1 1 , PM< l I , YM(  I )  I 

21  CONTINUE 

T  (  t  )*T<  D-TGII) 

0<  I  )>0< I )-QG( I ) 

NF<  I ) *NF I  I )-NFGI I ) 

SF< I  )  *SF(  I  )-SFG< I ) 

PM <  I  )«PM< | J-PMG ( I ) 

YM< I  |  *YM  ( I l-YMGI 1 ) 

20  CONTINUE 

IF  (NCASE.NE.C)  GO  TO  53 
DO  54  1*1. NPTS 

TGI  I ) *T(  I) 

QG( I )«0( I ) 

NFG(  I  )  *NF( I ) 

SFG(  I  )-SF< I ) 

PMG(  I  ) -PHI  I ) 

YMG<  I  )-YM< I ) 

54  CONTINUF 
NCASE  «  1 
53  WRITE  <6,821 

82  FORMAT  ( 27H1TARE  CORRECTED  PERFORMANCE) 

WRITF  <6,25) 

25  FORMAT  < IhO ,2X,  IHI . 6X. 5HTILB ) , 1  OX , BHQ IF  T-LB) »  7  X ,6HNF ( l B I ,9X  , 
16HSFILB),9X,9HPM<FT-LB) ,6X,9HYM(F T-LB) » 

WRITE  16,26)  (I,T<I),Q<l)«NFII),SF(l),PM<I),YM(I),I*ltNPTS) 

26  FORMAT  (I4.6E15.5/) 

CALL  HARMON  <T,18HPR0P  THRUST  LB  ,18,0, NPTS ,0 .0 ) 

CALL  HARMON  <Q,18HPR0P  TORQUE  FT-LB, IB ,t , NPTS ,0 . 0 ) 

CALL  HARMON  <NF,24HPR0P  NORMAL  FORCE  LB  ,24 ,0 ,NPTS , p  .0  ) 

CALL  HARMON  <SF,24HPR0P  SIDE  FORCE  LB  ,24,0, NPTS, 3.0) 

CALL  HARMON  <PM,30HPROP  PITCHING  MOMENT  FT-L8  ,30, C , NPTS, 3 .0  I 
CALL  HARMON  IYM.30HPROP  YAWING  MOMENT  FT-LB  ,30,0, NPT  3,0 .0 ) 

GO  TO  10C0 

2CG  RE AO  (5,981  ( <  SA  <  I  ,  J  I  ,  I -1 , 3)  ,  J*1 ,4 1 

98  FORMAT  (12F6.0) 

WRITE  (6,99)  << SA(  I, J), 1-1,31 ,J-1,4) 

99  FORMAT  <  17HOSA  <R 1 ,R 2 ,R3,R4 ) / 1X,4(3F6. 2 , 3X1 ) 

WRITE  (6,29) 

29  FORMAT  (3HCSO) 

DO  70  J-1,4 

DO  71  K*l,  3 

READ  (5,27)  <  SD<  I ,  J,K  )  , I -1 ,NPTS) 

71  CONTINUE 

WRITE  (6,61)  STRSS(J) ,( ( SO ( I ,J,K) ,1*1, NPTS) ,K*1 ,3) 

70  CONTINUE 

OO  30  J-1,4 

WRITE  16,33)  STRSS(J) 

33  FORMAT  (  1H1.A2/1X, IHI , 8X, 2HE 1 • 11 X.2HE2 , 1 1 X.2HE3 , 11 X.2HE2 , 1 IX, 2HEX , 
16X, 6HTHETAD, 5X,4HSIG2 ,9X,4HSIGX,5X,12HTS(LB/SQ  IN), IX, 

212HSSILB/SQ  INI) 
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00  31  1*1, NPTS 

00  32  K-1,3 

E(K)*SCC(K,JI*SA(K,  JI*(SD<  I,J,K)-SZD(  J,KII*ASIGN<  l,KI 
32  CONTINUE 

t>IEIll4EI3)l»0.5 

B*((E(  II-EI2I  )**2*(E(3)-E(2ll**2l*C.5 

C-SQRTIBI 

EZ«A*C 

EX«A-C 

ATN*2.0*IA-E(2I  » 

ATD*E( 3 1  —  E  < 1 ) 

THETA*ATAN2( ATN.ATOI 
THET AO* THETA* 5 7. 2957 79 5 
SIGZ*YMF*(  EZ*PR*EX I 
SIGX-YMF*(EX*PR*EZI 
SIGZPX*SIGZ*SIGX 
SIGZMX-SIGZ-SIGX 

TSI I,JI»(SIGZPX-SIGZMX*SIN( THE  TAM *C. 5 
SSU,J)«-  SIGZMX*COSITHETA)*0.5 

WRITE  (6,341  I,E(ll,F(?l»F(31*EZ,EX,THETAf),SlGZ*SIGX,TS(l*J), 

1SS( I,Jt 

34  FORMAT  (  /  IX,  1 2,  IX, 5E 13. 5,F8. 2 ,4E 1 3. 5  I 
31  CONTINUE 

30  CONTINUE 

IF  (NPTS.EQ.il  GO  TO  lClO 
00  35  J-1,4 

CALI  HARMON  (TSI 1 « J I « 24HTENSI LE  STRESS  LB/SQ  I  N,24, J, NPTS , RM( J II 

35  CONTINUE 

00  36  J*l,4 

CALL  HARMON  ( SS ( 1 ,  J I  , 30HSHE ARI NG  STRESS  LB/SQ  !N  , 33, J, NPTS, 
IRM(JI) 

36  CONTINUE 
GO  TO  10G0 
END 
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<*  o  o 


1BFTC  SSHARM  L  1ST, DECK, RFF 

PART  5  SUBROUTINE  HARMON 

HARMONICALLY  ANALYSES  A  FUNCTION  GIVEN  AT  IA  POINTS. MAX  IA-90 
SUBROUTINE  HARMON  I FN,BCD,NBCD,NS, I A,RM) 

DIMENSION  A I  44 1 ,B<44 I ,GAM( 44) , AAI44I  ,BCD( 12) ,CPHF (90,44) t 
l SPHE( 90*  44 ) ,FN( 9C  ) 

CATA  KOUNT/O/ 

IF  (KOUNT.F.O.l)  GO  TO  70 
PI  -  3.1415926536 
TP  I ■ 2.0*P I 
OPR  «  180.C/PI 
F  IA  -  FLOAT ( I  A ) 

TOA  •  2.0/FIA 
NFC  «  (  IA-D/2 
00  80  1-1,  IA 

F  I*FLOAT (  1 ) 

00  90  NN«l,NFC 
FNN«FLOAT(NN) 

XC«FNN*(FI-l.O)/F!A 
XC>AMOO( XC, l .0  ) 

PHE-XC*TPI 

CPHE ( 1,NN)-C0S(PHEI 

IF  (ABSICPHEI  I,NN)).LF.  .0001)  CPHEI I ,NN)«0.0 
SPHE ( I ,NN I «S! N( PHE I 

IF  (ABS(SPHE(  I.NND.LE.  .0001)  SPHE  (  I  ,NN)  >0.0 
90  CONTINUE 
8C  CONTINUE 
MOUNT  >  1 

70  00  10  NN* 1 ,NFC 
AINNI  >  0.0 
B(NN )  >  0.0 
00  2C  1-1,  IA 

A(NN)  -  A(NNt«-FN( I l*CPHE( I,NN) 

B(NN)  -  B(NN)+FN( I )*SPHE(  l,NN) 

20  CONTINUE 

AINN)  ■  T0A«A(NNI 
B(NN)  >  TDA*B( NN ) 

ARGAA  >  AINN )**2*B( NN 1**2 
AA(NN)  ■  SORT! ARGAA ) 

IF  (AINN) .NE .  0.0  .OR.  B(NN)  .NE.  0.0)  GO  TO  30 
GAM(NN)  ■  -0.0 
GO  TO  10 

30  GAM(NN)  -  - ATAN2 ( B ( NN ) ,A(NN) ) 

GAM(NN)  -  GAM(NNI*OPR 
10  CONTINUE 
AZERO  >  0.0 
00  40  1*1,  IA 

AZERO  *  AZEROtFN(I) 

40  CONTINUE 

AZERO  >  AZERO/F I A 
NBM-NBCD/6 

WRITE  (6,411  (BCO( I ), I-l.NBM) 

41  FORMAT! 1HI , 39X, 4CHF0UR I ER  COEFFICIENTS  AND  PHASE  ANGLES  OF, 

1 IX, 1 2A6)  ’ 


81 


IF  (NS.NE.O)  WRITE  (6,421  NS, RM 
42  FORMAT  I  18H  BLADE  SEGMENT  NO..I3/8H  RADIUS-, GI2. 5, 4H  FT . » 

WRITE  (6,511  AZERO 

51  FORMAT  ( //45X,8HA( ZERO  I *G 11,4/1 
WRITE  (6,501 

50  FORMAT  ( 1  3X, 14HHARM0NIC  ORDER, 10X ,4HA ( N) ,16X,4HB< N)  , 1 5X, 5HAA( N) , 
1 1 4X  *  6HGAM ( N ) / 1 
DO  60  NN-1,NFC 

WRITE  (6,521  NN,A(NN),B(NN) , AA (NNI ,GAM(NNI 
60  CONTINUE 

52  FORMAT  ( 18X, 12,4X,4( 1CX  ,G1 1 .4 > /> 

RETURN 
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